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A completely new mechanism to generate the observed amount of large-scale cosmological mag- 
netic fields is introduced in the context of three-form inflation. The amplification of the fields 
occurs via fourth order dynamics of the vector perturbations and avoids the backreaction problem 
that plagues most previously introduced mechanisms. 



PACS numbers: 98.80.Cq 

INTRODUCTION AND CONCLUSION 

Large scale magnetic fields exist in most observed 
structures in our Universe such as galaxies, clusters, fil- 
aments, and beyond, stretching to regions in the sky 
where little to next to no gas is present. Observational 
evidence is growing by the day, unveiling an unsettling 
uniform and ubiquitous presence; yet, their nature and 
origin remains draped by a hazy veil of mystery, see the 
reviews [IH1], and [SHS]- 

The property which is perhaps most intricate to frame 
within a consistent theoretical model for cosmic mag- 
netic fields is their pervasively wide coherence length, 
which can stretch well beyond the Mpc mark; cosmo- 
logical inflation might be the only viable candidate for 
its automatic transformation of short wavelength fluctu- 
ations into beyond-the-horizon correlations [ID] . 

Many simple options for inflationary magnetogenesis 
suffer from a backreaction issue, see [HUl4] . In short, as 
inflation proceeds at amplifying electromagnetic waves, 
the overall energy density of the latter catches up with 
that of the former, leading to instability in the de Sitter 
background. When this occurs inflation stops, and if this 
happens too early then we are not able to generate an 
almost scale invariant spectrum of curvature perturba- 
tions, as well as expand the causal region of the Universe 
beyond the observed one today - a possible viable model 
is perhaps the simplest [ID] , but it appears to be success- 
ful only in the strong coupling (hence non perturbatively 
trustable) regime [HI Q2] (but see the proposal of [T5] 
and criticism in [16]). 

The reasons behind this unhappy ending are twofold. 
In the first place, what is generally found is that the elec- 
tric field, while strongly suppressed in a highly conduc- 
tive plasma (such as the primordial soup during and after 
reheating) , plays the part of the giant in vacuum, see |13j . 
Secondly, the electric field tends to be redder than the 
magnetic field [13], and since inflation generally blindly 
processes an extremely wide range of modes, the mag- 
netic spectrum needs to be close to being red or all the 
allowed energy density will be snatched by small scales; 
but if this is so, the electric field will greatly dominate 



in the IR, leaving only crumbs of power to the magnetic 
field. 

In this work we show how these general issues can be 
circumvented in a particular model of inflation, three- 
form driven inflation [TTJ [18] , thanks to the peculiar re- 
shaping of the magnetic spectrum the coupling between 
the three-form and electromagnetism produces. 

Three-form inflation is a generally viable alternative 
to scalar field models. Details of the rich dynamics that 
three-form inflation allows can be found in various cases 
in the references [T7H22] ; for instance, the simplest mod- 
els cannot produce phantom inflation without instabili- 
ties or blue-tilted spectrum, but there are viable models 
generating slight red tilt with or without slow roll. 

Once coupled to EM, three-form inflation provides us 
with the two key ingredients for a promising magnetisa- 
tion of the large scale Universe. First, it entails the pos- 
sibility of exponential growth in the EM sector, driven by 
an instability in the three-form rotational sector. More- 
over, and more importantly, provides an effective UV cut- 
off at which the instability dies off; it then automatically 
splits it from the scale of inflation, allowing for a careful 
selection of the modes that are amplified, thereby upset- 
ting the democraticity of inflation. 

This can be understood in the dual description of 
three-form inflation as a scalar field inflation, since the 
duality maps kinetic into potential and vice versa. How- 
ever, the duality can break down in the presence of non- 
trivial self-interactions [22] of the three-form or nonmin- 
imal couplings to geometry [IT] or to matter [31]. Thus 
three-form inflation is equivalently described by a scalar 
field model only in the simplest minimally coupled cases. 
The novel aspect here is the coupling to the electromag- 
netic field: this both promotes the model to a physically 
distinct alternative to scalar-field driven constructions by 
breaking the duality (as it is apparent when one sees that 
the three-form has propagating vector degrees of freee- 
dom), and resolves a generic shortcoming of theirs by 
enabling sufficient magnetic field production without the 
backreaction issue. 

Our result are shown in Figs, [I] and [2] despite the 
electric field still being the one to watch out for, mag- 
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netic fields which today have strengths in the 10~ 15 
Gauss range can be easily produced without excessive 
fine-tuning of the parameters, and without jeopardising 
the development and completion of a long epoch of infla- 
tion. 



THE COUPLING 

Let be the photon vector potential and B pvp the 
three-form. The canonical lagrangian including the both 
fields is 

Ca + C b = -\f 2 (A) - lf 2 (B) - V(B 2 ) , (1) 

where the Faraday forms are computed from an n-form 
potential N as F(iV) Ml ... Mn+1 = (n + Fjld^iV^...^^]. 
The components of the dual of the three-form are [23] 



B a = ^af3 1 sB !3 ' lS . 



(2) 



The most general Lorentz-invariant, quadratic, second 
order and U{1) invariant coupling of the two fields is of 
the form 



Cab = 



-aF^A)F^(B) 



(3) 



THE AMPLIFICATION MECHANISM 



We decompose the spatial part of the vector potential 
A^ = (Ao,Ai) in terms of its transverse and longitudi- 
nal components as Ai — Af + di\, where V • A T = 0. 
The three-form is decomposed in a similar way where 
Bi — Bj + c*i£. Since the physical photons correspond 
to the transverse degrees of freedom, we are particularly 
interested in the vector perturbations. The vector part 
of the Einstein-Hilbert Lagrangian coupled to the three- 
form is EH 
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,(4) 



where we have used the background Friedmann equa- 
tion and C T is the gauge-invariant rotational degree of 
freedom of the metric perturbations [2 5) ; X is the back- 
ground value of the three-form field X = sf—B 2 in the 
notation of Ref. [52]. Eliminating the unphysical A 
through A n = %' — a(B — £') 26J, the relevant part 
of the action Q simplifies to 

1 
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(5) 



where i] is conformal time, and A = 8 lJ didj. The scalar 
polarisations of the photon are nondynamical, as ex- 
pected. It is interesting to note though that the scalar 



part of the interaction decouples from the electromag- 
netic field, it contributes to the effective sound speed of 
the three-form. 

The equations of motion for the vector degrees of free- 
dom become (omitting from now on the superscripts T, 
since all vectors are considered to be transverse) 

M|fc 2 C + 2V>(B-XC) = 0, (6) 
V x (B - XC) - aX(~dl + A) A = , (7) 



-dl + A)(A + 2aB) = 0. 



(8) 



We note that in the absence of the coupling, both the 
metric and three-form vector perturbations are nondy- 
namical. We can eliminate C and obtain a closed equa- 
tion for the three-form perturbation in Fourier space: 



where 



f = 2« 5 
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(9) 



(10) 



and B is the Fourier transform of B. Thus the three- 
form rotational modes propagate with a nontrivial dis- 
persion relation; they can therefore in principle be sig- 
nificant even at large scales. In particular, there is a 
divergence at / -> 0, that is, 2VxX -t Mpk 2 . Depend- 
ing on the sign of the sound speed squared, which in turn 
is determined by the shape of the potential, the vector 
perturbations over some range of scales can be drastically 
amplified (see below). 

Plugging the solution into equation ^ then shows that 
the electromagnetic potential evolves as a harmonic os- 
cillator driven by an external force, A" + k 2 A = F(B), 
where F(B) — —2a(d 2 + k 2 )B and A is the Fourier trans- 
form of A. One can also consider the autonomous evolu- 
tion equation for A, 



(-9 2 + A)[l-/(-^ + A)]A = 0. (11) 



From this it is obvious that the usual plane waves so- 
lutions always exist, but now there are two additional 
modes. At / ~ 0, where the three-form perturbation is 



divergent, the equation (111 reduces to second order. Go- 



ing to Fourier space, we have for each propagating degree 
of freedom 



A^ 



(/" - 1 + 2fc 2 /)^" 



2 ! -k 2 A' 



k\ 



f" - 1 
/ 



fc 2 )^ = 0. 



(12) 



Notice the appearance of fourth order time derivative 
in our equations, which might lead to the appearance 
of negative norm states. We are working here with an 
effective field theory which is only valid at low energies 
k <C Vq^ 4 , where Vq sets the UV mass scale of the model, 
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and we need to ensure that the dynamic effects are con- 
fined to such range; we will see that this is indeed the 
case in our specific examples. 

Hence, to summarise, the effective driving force for A 
induced by the instability of B provides a mechanism to 
boost the magnetic fields. We are not aware of any pre- 
vious amplification mechanism based on nontrivial dy- 
namics of the primordial rotational perturbations, since 
most magnetogeneses in the literature can be effectively 
described by introducing either a mass term or a time 
dependent coupling constant to the photon [llj . 

EXAMPLES 

For simplicity, we consider de Sitter solutions with con- 
stant comoving field X. There are two classes of such 
fixed points [55]: class B, which are the critical points 
corresponding to X 2 = ±|Mp, and class C, which are 
the critical points corresponding to the minima of the 
potential. The stability of these fixed points depends on 
the shape of the potential. In both cases the evolution 
equation dl2|) simplifies to 



lr { ^-k 2 )A = 0, (13) 
Jo 



- f I - 2k 2 ) A" 



f is given by Eq. (10 1 evaluated at the fixed point. We 



can readily solve this equation to obtain 

A(t}) = Ai cos (kr)) + A 2 sin {kr}) + A 3 e rkv + ^ 4 e _r fl'4) 

where Aj are constant vectors and T 2 + 1 = l/(fok 2 ). 
Besides the usual plane waves, the two interesting and 
new solutions are also oscillatory if T 2 < 0, but if T 2 > 
0, one of them is exponentially growing and the other 
exponentially decaying. 

As an example consider the exponential model V — 
V exp (— f3X 2 /Mp) . At the fixed point B, which is stable 
if /3 is positive, we have 



(15) 



A 2 4 



where we have defined A 2 = 8a 2 /3/(l - 8a 2 /3) ~ 8a 2 /3 
and k 2 A = 8/3U/(3M 2 A 2 ) ~ /3V7(a 2 M 2 ). In this way it is 
immediate to see that the exponentially growing solution 
exists only for k < k\; k\ is therefore the UV scale which 
defines instability band. We will repeatedly be using the 
notation k = k/T-Lf, which normalises k to the highest 
mode accessible by inflation, T-L e - we denote T-L = a' /a, 
where a is the scale factor of the Universe and the prime 
stands for conformal time derivative. 

Therefore, the scale acting as de facto UV cutoff is 
identified with, &a; this is well below the actual UV 
cutoff, set by V^ 4 , and our theory is trustable within 
this limit. In the more general and complete theory, 
alongside terms in the Lagrangian like F(A)F(B), which 



would make the Hamiltonian unbounded from below, also 
higher order, stabilising, terms such as, schematically, 
(F(A)F(B)) 2 /Vo would appear; in the full theory, if con- 
structed appropriately, there will be no ghosts (see for 
instance the Galileon case) - we believe that it goes be- 
yond the scope of the present work to show explicitly that 
such a theory can be assembled, and we shall be content 
with working in the, safe, low energy limit. 

The other fixed point C at the origin is an attractor 
when ft is negative. There we have 



2;„2 



T z k 



a 2 M 2 



(16) 



Thus, when this point is an attractor, there is no insta- 
bility of the vector modes. If we are inflating at a saddle 
point (i.e., j3 > 0) vector modes will be unstable at scales 

We summarise the properties of the fixed points for 
several classes of potentials and the corresponding ex- 
pressions for r 2 + 1 in table [i] 

In general, it is easy to see when the instability for the 
vector potential A, parametrised by T 2 , is present, and 
at which scales. If large scales, k — > 0, are to be exponen- 
tially amplified, then it is necessary that 4a 2 A 2 /Mp < 1. 
Notice that at large momenta r 2 — > — 1 and we recover 
plane wave solutions. Hence, there is in general a low en- 
ergy band which is unstable given the above condition. 

Beside the turning point where T 2 changes sign, typi- 
cally there would be a singularity in T 2 for the mode k 
which nullifies its denominator. This can lie within the 
exponentially growing band resulting in an infinite en- 
ergy density around such scale. The question of where 
the singularity will occur for a specific potential depends 
on the sign of VxX: if this is positive the divergence will 
appear; in the opposite case VxX < the problematic 
scale is beyond the UV cutoff. In the specific examples we 
discuss this uncontrolled instability is avoided, whereas 
the exponential, controlled, growth is attained for a low 
energy band of modes. 



THE SPECTRUM 

The action ^ is canonical for the EM field, and as long 
as we treat the coupling to the three form as a perturba- 
tion we can quantise canonically [57]. We are interested 
of course in solutions which grow with conformal time, 
so we will choose A\ = A% = A4 = - although this 
may appear artificial, it is immediate to convince one- 
self that, were these solutions kept, they would become 
very rapidly irrelevant in comparison to the growing A3 
term (Ai and Ai are simply vacuum plane waves, and 
A4 decays exponentially). 

We have two options in fixing the initial conditions (IC) 
for the EM vector potential A in Fourier space: either we 
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TABLE I. Stability of the de Sitter fixed points and the characteristic width of instability of the vector potential in four classes 
of models. U - unstable; S - stable. 



set them at some initial time, which can be taken to cor- 
respond to the beginning of inflation r]i, for all modes, 
or we can normalise the fields to their Bunch-Davies vac- 
uum V2kA = eT lkri at the time they become superhori- 
zon. The first possibility is the simplest, but perhaps not 
the most physical nor appealing one, for it introduces an 
explicit and very strong dependence on the IC of infla- 
tion, which is antithetic to what inflation is brought in 
for. Also, it poses the theoretical question on how to 
properly choose a vacuum state since the effective cou- 
pling is never turned off. Nonetheless, the latter issue 
can be resolved if we work in the small coupling regime: 
in this case it is consistent to choose the Bunch-Davies 
solution as the approximate IC. 

The other choice corresponds to a coupling which dy- 
namically becomes relevant at some time (e.g. horizon 
crossing), while being inefficient for smaller scales. This 
option makes it clear that for all modes the specific IC do 
not matter, while they can be consistently chosen to be 
in vacuum at any early time without affecting the final 
result. However, this demands a more elaborate theoret- 
ical motivation for the time-dependence in a - if possible 
at all - and will be discussed in a subsequent publica- 
tion gg. 

One more comment is in order here. Since we work 
with the simplified and idealised toy model where X does 
not go anywhere, also / is constant, and the equation of 
motion for B ^ does not admit plane wave solutions at 
any time. However, it is easy to convince ourselves that 
in more realistic models, when fk 2 — > ±oo as krj — > — oo 
(which depends on the shape of the potential) then the 
rotational perturbation B is initially in its Bunch-Davies 
vacuum, and so will A. That is, the exponent T in re- 
ality is time-dependent as well, and flips from imaginary 
to real at a given time; we take the latter to be rji for 
simplicity. This also means that in the full model the de- 
pendence on the initial conditions will be pronouncedly 
less severe. 

The IC fix the value of A3 and give us the following 
solution to work with: 



We can immediately write down the corresponding ex- 
pressions for the energy density associated with the EM 
field at the end of inflation as 



Pem = 



4tt 2 

nji 

8ir 2 at 



(18) 

dn K 3 (r 2 + 1) e 2r «(V«*-i) , (19) 



where of course n e = 1 by definition, and we can clearly 
ignore the second piece in the exponent, being much 
smaller than the first, positive, contribution. Also, we 
work with small a, which means Ac 1 and r> 1; this 
implies that in general E 2 3> B 2 . 



The integral (191 is not simply dealt with, and we must 
resort to numeric methods to solve it. The limits of in- 
tegration are n = K{ = exp(— N), where N is the total 
number of e-folds inflation lasts for, and n = k\, beyond 



A 



1 



ji-fcOh-i?) IC at ry = ?7i . 



(17) 



which the solution (17) does not apply anymore. 

In Fig. [l] we present the regions in the parameter 
space (A, k^) for the well-behaved exponential potential 
exp(— /3X 2 ), for which the ratio pem/px stays below one 
at the end of inflation (blue region) and goes beyond one, 
thereby leading to instabilty in the de Sitter background 
(red region). Superimposed are the contours for a given 
value of the magnetic power at 3/Gpc today, where the 
regions for which 5^ > 10 -15 Gauss and Sg > 1CT 25 
Gauss are in evidence (green and yellow, respectively) 
- we define (Sg) 2 = k 5 \A\ 2 / 4ir 2 . These are values that 
correspond to the magnetic fields observed in the inter- 
galactic medium, and a typical value for a successful seed 
to be fed the magnetohydrodynamic plasma at late times. 

Fig. ^ plots the power spectrum today, for A = 10~ 2 
and fcA/fcmin — 10 2 ; notice the knee at high energy and 
the following very rapid decay. Thanks to the UV cut- 
off, such a spectrum also nicely helps in avoiding the 
stringent constraints which late time cosmology poses on 
smaller scale fields [55]. These figures were obtained for 
around N — 66 e-foldings of inflation, H e /a e ~ 10 13 
GeV, and assuming the Universe was always dominated 
by radiation from the end of inflation onwards. 

There is some degree of fine-tuning: as one realises 
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FIG. 1. Allowed (in light blue) and disallowed (in light red) 
regions in the (log A, log (fcA/fcmm)) parameters space. Super- 
imposed are contours of constant log S% today, at 3/Gpc (in 
Gauss units). 
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FIG. 2. The power spectrum today, 6% (Gauss), for A = 10 2 
and fcA/fcmin = 10 2 . 



by looking at Eq. (17 1, these is a very strong sensitiv- 



ity (doubly exponentially in terms of the total number 
of e-folds N) to the initial time, and only scales not too 
far from the first ones available to inflation can be am- 
plified efficiently without leading to a catastrofic back- 
reaction effect. In the figure we show the final power 
at k — 3/Gpc, where the largest mode is only around 
1 /Hq « 4Gpc. The plot of the spectrum makes it clear 
that only a couple of orders of magnitude in k will be am- 
plified efficiently; for instance, the strength at 10/Gpc is 
10~ 21 Gauss, but drops to 10 -47 Gauss at just 30/Gpc. 
If we were to push the beginning of inflation further back, 
at much smaller Hi, it would be impossible to amplify, 



without backreacting, today's large scales. With the al- 
ternative prescription of a time- dependent coupling the 
sensitivity to initial conditions can in principle be elim- 
inated. Recall again also that this feature is partly due 
to the idealisation of the model we work with. 

Looking again at the allowed region in Fig. [T] it ap- 
pears that one may push the coupling a ~ A to arbitrar- 
ily small values and still come out with enough power 
on a given scale. This is only partially true, since a 
also enters the UV cutoff, and as a — > one needs to 
keep a 2 k\ ~ /3Poe~ 2/3 / 3 /M| fixed, that is, the window 
of modes which can be boosted has to be kept narrow. 
Order one slopes /3 would then demand extremely low 
bare Vo, which is desirable since it corresponds to little 
kinetic energy in the dual scalar description. Further- 
more, one can consider models with either very high or 
very low potential slope parameters f3 to reconcile suit- 
able UV cutoffs with small couplings a. 



ADDENDUM 

After this work was completed, two preprints dealing 
with the anisotropics generated in inflationary magnc- 
togenesis appeared 29, 30 . Their conclusion is that 
such contributions to curvature perturbations become 
unacceptably large during the subsequent radiation era, 
thereby lapidarily coercing largest classes of models to 
failure. Although we have not completed the full, and 
very laborious and lengthy calculation, we believe that 
this result needs not apply as is to our mechanism for 
three reasons. 

The result of [29l [30] strongly depends on the limita- 
tions on the generated power spectrum; our case presents 
a very peculiar spectrum, due to its low energy cutoff, 
which allows for almost any slope, thereby offering the 
possibility that the dangerous term in fact be not as big 
as thought. 

Secondly, our Bardeen equation is quite different due 
to the three- form perturbations [22] ; the solutions to this 
equation during inflation can have different scaling with 
momentum k and defuse the too large anisotropies. 

Finally, the time evolution of A is not a power law 
as in [211 ED], but exponential, which in turn is strongly 
fc-dependent; this implies that the time-evolution of per- 
turbations could further differ from what predicted in 
their case. 

We elaborate on this critical issue in our upcoming 
publication [23]. 
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